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00 Abstract 

We study a supersymmetric version of the seesaw mechanism type-Ill considering two variants 

Qh of the model: a minimal version for explaining neutrino data with only two copies of 24 superfields 

Q_i and a model with three generations of 24-plets. The latter predicts in general rates for u — > ej 

inconsistent with experimental data. However, this bound can be evaded if certain special conditions 

within the neutrino sector are fulfilled. In case of two 24-plets lepton flavour violation constraints 

can be satisfied much easier. After specifying the corresponding regions in the mSugra parameter 

ON space we show that under favorable conditions one can test the corresponding flavour structures in 

OO 

£\j the leptonic sector at the LHC. For this we perform Monte Carlo studies for the signals taking also 

into account the SUSY background. We find that it is only of minor importance for the scenarios 
studied here. 
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I. INTRODUCTION 



Neutrino oscillation experiments give currently the main indication for physics beyond 
the Standard Model (SM). The observed tiny neutrino masses can be easily explained by 
the seesaw mechanism, which at tree level can be written in just three different variants [lj, 
classified according to the SU(2) x £7(1) representation of the postulated heavy particles: 
type I postulates fermionic gauge singlets [2HH], type II scalar 5'[/(2)-triplets with hyper- 
charge 1 O [7] and type III fermionic triplets in the adjoint representation of SU(2) [8j. At 
low energies they all lead to a unique dimension-5 operator |9j HD] 

{m»U= f -f{HL){HL). (1) 

Neutrino experiments determine only f a p / A, but contain no information about the origin of 
this operator, nor about the absolute size of A. If / is a coefficient 0(1), current neutrino 
data indicates A < O(10 15 ) GeV. This value is close to, but slightly below the scale of 
hypothetical great unified theory (GUT) which should be larger than roughly 10 16 GeV to 
avoid bounds from the non-observation of proton decay. 

A possible way to stabilise the large hierarchy between the GUT scale and the electroweak 
scale is supersymmetry In its minimal form it leads to an unification of the gauge 

couplings in contrast to the SM [ T2HT8] . Moreover it can explain electroweak symmetry 
breaking as a radiative effect |19| . Supersymmetric variants of the different seesaw models 
have been considered by several authors, see e.g. |2UH25j . In these models renormalisation 
group evolution induce non-zero flavour mixing elements in the mass parameters of the 
sleptons even if they are flavour diagonal at the GUT scale. These in turn leads to sizeable 
contributions to lepton flavour violating (LFV) observables [26J. In case of seesaw type-I, 
low energy LFV decays such as U — > lj + 7 and Zj — > 31 j have been calculated in [2TH36] ; 
— e conversion in nuclei has been studied in [371 138] . 

To maintain gauge coupling unification the seesaw particles need to be included in com- 
plete SU(5) representations, i.e. one needs a 15-plet in case of type II and at least two 
24-plets in case of type III models. If one were to use only one 24-plet, then one would 
need either non-renormalizable operators at the GUT-scale (39] or an extended SU(5) Higgs 
sector (10] to explain neutrino data. The type-II and type-Ill models have received less at- 
tention than type-I. Note, however, that the former has actually fewer free parameters than 
type-I implying that ratios of LFV decays of leptons can actually be predicted as a function 
of neutrino angles in minimal supergravity (mSUGRA) as discussed in [21, 23J. For type-Ill 
it has been shown in ref. (25] that a generic model with three 24-plets is heavily constrained 
by the bounds on rare lepton decays, in particular due to the stringent bound on fi — > cy. 
The impact on — > in a gauge mediated supersymmetry breaking (GMSB) embedding 
of type-Ill was studied in ref. (41] . while in ref. (42] also possible LHC phenomenology from 
lepton-flavor violation has been discussed for the mSUGRA case. 

In this paper we are first going to show under which conditions the type-Ill model is con- 
sistent with the experimental data. This will then be compared with a two generation model 
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where the bounds due to /i — > are less stringent. Finally we will address the question 
to which extent LHC might observe lepton flavour violating processes in supersymmetric 
(SUSY) cascade decays. Compared to the previous study in ref. |42| we do not only consider 
the case of two as well as of three generations of 24-plets, but we also take into account the 
recently measured reactor angle #13. Moreover, we demonstrate in a Monte Carlo study for 
the LHC signal that the SUSY background is well under control. 

For the particle content we will assume the MSSM as framework. The recent observation 
of a Higgs like state at the LHC with a mass around 125 GeV |4"3"tl4*4"] can hardly be explained 
within GUT models with universal boundary conditions |4"5ll4"6]. The same holds in variants 
including seesaw states at high scales [47J. However, it is well-known that singlet extension 
of the MSSM can more easily explain a Higgs mass of this size as there are additional F- 
term contributions already at tree-level, see e.g. |4"8H50] and references therein. Including an 
additional singlet does not lead to any significant changes in the slepton sector as the singlet 
Yukawa couplings enter only at two-loop level in the RGEs of the corresponding parameters. 
For this reason we do consider bounds from direct searches at the LHC but do not take into 
account the requirement of correctly explaining a Higgs mass of 125 GeV. To explain also 
this condition and taking into account the theoretical as well experimental uncertainty it 
would be sufficient to shift the tree-level Higgs mass for all points in the following by about 
5 GeV. In the next-to minimal supersymmetric standard model (NMSSM) the tree-level 
mass of the light Higgs is given by |5T] . 

m% (cos 2 2(3 + ^sin 2 2p\ (2) 

where we introduced the superpotential coupling A of the singlet to the Higgs doublets. 
Thus, assuming tan = 10 one would need A ~ 0.6, i.e. not too close to the perturbativity 
bound of 0.75. However, the exact value of tan/3 plays only a sub-dominant role for our 
analysis in the following. If we choose tan /3 = 5, already A = 0.45 would be sufficient. 
Since the singlino couples only very weakly to the sleptons its role is negligible as long as it 
isn't the LSP. However, it can only be the LSP if the trilinear self coupling k of the singlet 
field is much smaller than A |52j . For example in the scale invariant NMSSM a bino LSP 
is a common feature j53]. Moreover, this constraint can easily be satisfied in more general 
singlet extensions with an explicit bilinear singlet term in the superpotential [541 [55] . 

This paper is organised as follows: in the next section we summarise the main features of 
the two variants of the type-Ill model. In section HI we first discuss how to accommodate 
the rare lepton decays in type III seesaw models. Afterwards we discuss lepton flavour 
violating signals from SUSY cascade decays and present the results of a Monte Carlo study. 
Finally we draw in section IV our conclusions. 



II. MODELS AND SPECTRA 

In this section we briefly summarise the main features of the supersymmetric version of 
the seesaw type-Ill model. In order to maintain gauge coupling unification for type-Ill we 



3 



add at the seesaw scale(s) additional particles to obtain complete S77(5) representation, i.e. 
a 24-plet. Note that the 24-plet actually also includes a gauge singlet and, thus, one has 
always a combination of the type I and the type III seesaw in this model. 

In the subsequent sections we present the superpotentials and the relation of the pa- 
rameters to neutrino physics. In addition, there are the corresponding soft SUSY breaking 
terms which, however, reduce at the electro-weak scale to the MSSM ones and, thus, are 
not discussed further. There are additional terms of the soft SUSY breaking potential, due 
to the heavy particles, that we do not discuss either, as their effect is at most of the order 
Mewsb/ M seesaw and, thus, can be safely neglected. 

In this paper we will assume common soft SUSY breaking parameters at the GUT-scale 
Mqut to specify the spectrum at the electro-weak scale: a common universal gaugino mass 
M1/2, a common scalar mass m and the trilinear coupling A Q which gets multiplied by the 
corresponding Yukawa couplings to obtain the trilinear couplings in the soft SUSY breaking 
Lagrangian. In addition the sign of the \x parameter is fixed, as is tan/3 = v u /vd (at the 
electro- weak scale), where Vd and v u are the the vacuum expectation values (vevs) of the 
neutral component of Hd and H u , respectively. The models discussed below also contain 
new bilinear parameters in the superpotential leading to additional bilinear terms in the 
soft SUSY breaking potential which are proportional to B Q of the MSSM Higgs sector. The 
corresponding RGEs decouple and their only effect is a small mass splitting between the new 
heavy scalar particles from the new heavy fermionic states of the order Bq/ M seesaw . This 
leads to a tiny effect in the calculation of the thresholds at the seesaw scale(s) [56] which, 
however, we can safely neglect. 

A. Supersymmetric seesaw type-Ill 

In the case of a seesaw model type-Ill one needs new fermions £ at the high scale belonging 
to the adjoint representation of SU(2). This has to be embedded in a 24-plet to obtain a 
complete 577(5) representation. The superpotential of the unbroken £77(5) relevant for our 
discussion is 

W = V2 5 M Y 5 W M 5 H - ^10 M Y 10 10 M 5tf + 5 H 2A M Y^ II 5 M + ^24 M M 24 24 M • (3) 

Here we have not specified the Higgs sector responsible for the S77(5) breaking as it only 
enters logarithmically via threshold corrections at the GUT-scale and, thus, plays a minor 
role for the subsequent discussion. The new parts, which will give the seesaw mechanism, 
come from the 24m- It decomposes under S77(3) x SU(2) x U(l) as 

24 A/ = (1, 1, 0) + (8, 1, 0) + (1, 3, 0) + (3, 2, -5/6) + (3*, 2, 5/6) , (4) 

= B M + Gm + Wm + X M + X M ■ 

The fermionic components of (1,1,0) and (1,3,0) have exactly the same quantum numbers 
as a right-handed neutrino u c and the required S77(2)-triplet E. Thus, the 24 M always 
produces a combination of the type-I and type-Ill seesaws. 
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In the SU(5) broken phase the superpotential becomes 

Wui = W M ssm + H U (W M Y N - ^B M Y B )L + H U X M Y X D C 

+^B M M B B M + -d M M G d M + ^WmM w Wm + X M M X X M .s (5) 

As before we use at the GUT scale the boundary condition Yn = Yb = Yx and Mb = 
Mq = My/ = Mx- Ypf, Yb and Yx are n x 3 while Mq, My/ and Mx are n x n-dimensional 
matrices if we include n generations of 24-plets. Integrating out the heavy fields yields the 
following formula for the neutrino masses at the low scale: 

m » = ~f {™*» M ? Yb + l Y w M wYw^j ■ (6) 

As mentioned above there are two contributions, one from the gauge singlet the other from 
the SU(2) triplet. In this case the calculation of the Yukawa couplings in terms of a given 
high scale spectrum is more complicated than in the other two types of seesaw models. 
However, as we start from universal couplings and masses at Mqut we find that at the 
seesaw scale one still has Mb — My/ and Yb — Yy/ so that one can write in a good 
approximation 

m u = -vl^Y?M w l Y w . (7) 

Being complex symmetric, the light Majorana neutrino mass matrix in eq. ([7]), is diago- 
nalized by a unitary 3x3 matrix U [B] 

m v = U T -m v -U . (8) 

Inverting the seesaw equation, eq. (|7j), allows to express Y w as [57] 

Yy/= % ^sJW M -R-^a v -V\ (9) 

for n = 3 where the fa v and Wm are diagonal matrices containing the corresponding eigen- 
values. R is in general a complex orthogonal matrix which is characterised by three angles 
<pi which are in general complex. Note that, in the special case R — 1, Yy/ contains only 
"diagonal" products \J M^rrii . For U we will use the standard form 

' C12C13 S12C13 s l3 e~ tS \ 

U = 



-S12C23 - Ci 2 S23Sl 3 e* 5 C12C23 - Si2S 2 3Sl 3 e i ' 5 S 23 Ci 3 
\ S12S23 — Ci2C23Si3e 4<5 — Ci2S2 3 — Si2C23Si 3 e J<5 C2 3 Ci3 / 



/ e iai/2 q \ 

e ia2 ' 2 

V V 



(10) 



with Cij = cos(9ij) and Sy = sin(^y). The angles 612, ^13 and 623 are the solar neutrino angle, 
the reactor (or CHOOZ) angle and the atmospheric neutrino mixing angle, respectively. 5 
is the Dirac phase and a, are Majorana phases. In the following we will set the latter to 
and consider for 5 mainly the cases and 7r. 
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B. Supersymmetric seesaw type-Ill with two 24-plets 



Current neutrino experiments only determine the differences of the neutrino masses 
squared. Thus it might well be that only two of the light neutrinos are massive whereas 
the third is either massless or has a mass much smaller than the others. Such a situation 
is obtained if only two 24-plets are present similarly to the case of the seesaw type I with 
two right-handed neutrinos as discussed for example in |5TM)U] . We call this class of models 
3x2 seesaw, see also |61| . 

In the following we work in the basis where Mw is a 2 x 2 diagonal matrix denoting the 
eigenvalues by M{ (i = 1,2). Similarly to the three generation case one can express the Yw 
in terms of low-energy neutrino parameter and model dependent high energy parameters as 
also discussed in the context of seesaw I models |62|: 



Y w = \l-^M^R(VfK l ]u\ (11) 



2 v u 

The i?-matrix is now 2x3 which can assume the following forms: for normal hierarchy 
(mi = 0) 

_ fo cos(0) -sin(0)\ 
n norm - y Q g . n( ^ cog( ^ j [LZ) 

in case of normal hierarchy in the neutrino sector (mi = 0) or and for inverse hierarchy 
(m 3 = 0) 

_ / cos(0) - sin(0) ()\ 
norm ~ \sm((f>) cos(0) J [ } 

in case of inverse hierarchy (777.3 = 0). Note that the i?-matrix is parametrised by one 
complex angle only in contrast to the three generation case. 



C. Lepton flavour violation in the slepton sector 

From a one-step integration of the RGEs one gets, assuming mSUGRA boundary con- 
ditions, a first rough estimate for the lepton flavour violating entries in the slepton mass 
parameters: 

(Aml)^-^- 2 (3ml + Al)(Y^LY^) , (14) 
{AA) ltij ~ ~a k -^A (Y e Y#LY*) , (15) 

l07T z V / ij 

for % 7^ j in the basis where Y e is diagonal, = \tl(Mgut / 'Mi)8ij and Y^ is the additional 
Yukawa coupling where k indicates the number of 24-plets. 

6 9 

a-i — - and a 3 = - . (16) 
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Both models have in common that they predict negligible flavour violation for the right- 
sleptons 

(Am|) y ~ (17) 

which is a general feature of the usual seesaw models |2o] . Although it is known that approx- 
imations eqs. (14) and (15) do not reproduce well the actual size of the off-diagonal elements 
they do give the functional dependencies on the high scale parameters. Therefore they are 
a useful indicator on how the rare lepton decays U — > Ijj depends on these parameters as 
the corresponding decay modes scale roughly like 



| (Am 



12 

Br{k Ijj) oc a 3 mf. IV "Z^ tan 2 /3. (18) 



rn 



where fh is the average of the SUSY masses involved in the loops. Using the parametrization 
for the Yukawa couplings of eq. (pi the entries in (Am^jy can be expressed as 



(AmDij oc UiJJjpy/fn^y/fn^I^RkpMk log (J^J ■ ( 19 ) 
In the special case that the matrix R is the identity matrix, eq. (19) reduces to 

(A?77|) 12 (X C 12 C 13 (~S 12 C 23 - Ci2S23Sl 3 e" J<5 ) Z\ (20) 
+ S12C13 (C12C23 - S 12 S 23 Si 3 e~ %& ) z 2 + S 23 C 13 S r3 e~ lS Z 3 

(Am 2 L ) 13 oc C12C13 (si 2 s 2 3 - ci2C2 3 si 3 e" i5 ) z x 

+ s 12Cl3 ( — C12S23 — Si2C23Si3e lS ) Z 2 + C2 3 Ci3Si3e lS Z3 
(Am 2 L ) 23 OC (S12S23 - Ci 2 C2 3 Si3e^ 5 ) (-si 2 c 23 - Ci 2 S23Sl3e J<5 ) z x 

+ (-C12S23 - s 12 c 23 s 13 e~ tS ) (c 12 c 23 - si 2 s 23 s 13 e tS ) z 2 



+ S23C23Ci 3 ^3 



where 



^log(^). (21) 



For the ansatz of degenerate seesaw states the combination M i \og(^ L ) becomes an overall 
factor, i.e. for degenerate Mb = M\y one may simply make the replacement Z{ — > rrii in 
eq. (20). For strict normal hierarchy, the expressions become even simpler. For instance, 
Am 2 L becomes 

(Am|) 12 oc (s 13 S23v / A(mi tm ) + A(m|) - ^ A{ml)s\ 2 + c l2 c 23 e iS ^ A{m%)s l2 ^ , (22) 

Inserting the best fit point data for oscillation parameters, except for S\s, and assuming 5 = n 
one can calculate the value for s 2 3 for which Am 2 L approximately vanishes as sf 3 = 0.0077, 
which agrees very well with the full numerical calculation shown in the next section, see 
Fig. U 
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Similar analytical estimates can be calculated in other limits and, even though absolute 
values for LFV processes are only rough estimates, ratios of LVF quantities can be calculated 
quite accurately in this way. 

In the numerical studies we will use the complete formulas as given in [52"| 155] . We will 
also consider the three body decays BR(7j — > 31 j) where we use the formulas given in |32] . 



III. NUMERICAL RESULTS 



In this section we present our numerical calculations. All results presented below have 
been obtained with the lepton flavour violating version of the program package SPheno 
[6^, 165] . The RGEs of the two seesaw III models have been calculated with SARAH [66-69J. 
For the Monte Carlo studies below we have used the SUSY Toolbox [70] to generate the 
interface to WH IZARD [7J]. 

All seesaw parameters as well as the soft SUSY breaking parameters are defined at Mqut- 
We evolve the RGEs to the scales corresponding to the GUT scale values of the masses of 
the heavy particles. The RGE evolution implies also a splitting of the heavy masses up to 
20% between the gauge singlet and the color octet. We therefore add at the corresponding 
scale the threshold effects due to the heavy particles to account for the different masses as 
discussed in [23]. However, since the gauge singlet doesn't contribute to the running of the 
gauge couplings, the main impact on gauge coupling unification is due to mass splitting 
between the color octet and the SU(2)l triplet. This splitting is for a seesaw scale of 
O(10 14 GeV) of the order of 10% and would result in a marginal shift of O(10~ 4 ) for the 
gauge couplings. Off-diagonal elements are induced in the mass matrices of the 24-plets. 
This implies that one has to go the corresponding mass eigenbasis before calculating the 
threshold effects. We use two-loop RGEs everywhere except stated otherwise. 

Unless mentioned otherwise, we fit neutrino mass squared differences to their best fit 
values [72] . Our numerical procedure is as follows: inverting the seesaw equation, see eqs. ^ 
and (11), one can get a first guess for the Yukawa couplings for any fixed values of the light 
neutrino masses (and angles) as a function of the corresponding triplet mass for any fixed 
value of the couplings. This first guess will not give the correct Yukawa couplings, since the 
neutrino masses and mixing angles are measured at low energy, whereas for the calculation 
of m v we need to insert the parameters at the high energy scale. However, it can be used to 
run numerically the RGEs to obtain the exact neutrino masses and angles (at low energies) 
for these input parameters. The difference between the results obtained numerically and the 
input numbers can then be minimized in a simple iterative procedure until convergence is 
achieved. As long as neutrino Yukawas are |l"V,y| < 1 Vi, j we reach convergence in a few 
steps. 
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FIG. 1. BR(^x — > cy) versus the seesaw scale for the two (dashed) and the three (solid) 24-plet 



scenario; mSUGRA parameters like in eq. (23), neutrino data fixed to the current experimental 
values, including sin 2 #13 = 0.026. 



A. Bounds from lepton decays 

It is known for some time that generically the supersymmetric seesaw III model predicts 
rates for fi — > e-y which are too large [25] to be compatible with the experimental bound 
for BR(,u — > ery) < 2.4 • 10 -12 |73| . However, this does not completely exclude the model as 
there are certain parameter regions where cancellations between different contributions can 
occur. In this section we explore the different possibilities. For the corresponding regions the 
question arises if they can be probed by other experiments, in particular the LHC From the 



discussion in the previous section, in particular eqs. (14) and (18), the rare leptons decays 
are mainly governed by the overall SUSY mass scale and the lepton flavour entries. The 
LFV entries in the softs are nearly completely governed by the choice of parameter in the 
heavy seesaw sector, while the dependence on the soft SUSY parameters is much weaker. 
We therefore fix the later to: 



m = 1000 GeV , M 1/2 = 1000 GeV , A = and // > . 



(23) 



In fig. [T] we recall the generic situation for the type-Ill seesaw. The dashed and full 
line correspond to the 2- and 3-generation model, respectively. Only a certain range for 
^Seesaw is allowed. The lower bounds stem from the fact that the gauge couplings become 
non-perturbative at the GUT-scale whereas the upper bounds are due to non-perturbative 
Yukawa couplings at the GUT-scale [25J. Every 24-plet contributes with Abi = 5 to the 
beta functions of the gauge couplings and, thus, obviously the possible range is larger for 
the 2-generation case compared to the 3-generation case. 

Equations ^ and (11) imply that one can induce special features for the Yukawa cou- 
plings when varying sin# 13 , the CP-phases and/or elements of the i?-matrix as has also been 
noted in ref. [31] in case of supersymmetric seesaw I models. This has an immediate impact 
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FIG. 2. BR(^x — > cy) (blue), BR(t — > ej) (red) and BR(r —> ^) (yellow) over the reactor angle 
013 for real parameters with Dirac phases 5 = (upper left), 5 = n (upper right) and 5 = 37r/4 
(lower panel). We set Mw = 10 14 • II3 and the SUSY parameters as in eq. (23). The dashed line 
indicates the current best-fit value for #13. 
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on the flavour mixing entries of the slepton mass parameters as can be seen from eqs. (14) 
and (15). As an example we show in fig. [2] the dependence on # 13 in the range allowed 
before the results of DAYA-BAY [74J and RENO [73] assuming three different values for 
the Dirac-phase 5 and a degenerate mass of 10 14 GeV for the 24-plets. Note that in this 
particular case the elements of i?-matrix do not play any role. As can be seen, 5 has to be 
close to 7r in this case to get BR(/i — > e^y) below the experimental bound. For completeness 
we note that the small spikes in the plots are numerical artifacts of our iterations procedure. 

With the recent measurement of #13 by Daya Bay and RENO one can now relate seesaw 
parameters from the requirement to respect the bound on /i — > cy^ As an example we fix 
in fig. [3] 5 = 7r, Mi = M 2 = 10 14 GeV and take R = II3. In this case the bound on /i — >■ ej is 
satisfied if M3 is close to 5 • 10 13 GeV. Note however, that the numbers obtained depend on 



1 For a recent update on /.t 
see [75]. 



ej in seesaw type I models taking into account the measured value of #13 
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FIG. 3. BR(/i -)• e7) (blue), BR(r -)• e7) (red) and BR(r -)• /ry) (yellow) as function of M 3 with 
#13 at the current best-fit value. We have taken 5 = 7T, Mi = M2 = 10 14 GeV and the other 



parameters as in eq. ( 23 ) 



PQ 




sin (0 13) 

FIG. 4. BR(/i — >■ e7) (blue), BR(t — > ej) (red) and BR(r — > ^7) (yellow) as function of sin 2 #i 3 
in the 2-generation model for 5 = vr, Mi = M 2 = 10 14 GeV (solid lines) and Mi = 2 • 10 14 GeV, 
10 14 GeV (dashed lines). The other parameters are as in eq. (23). 
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the SUSY point chosen in parameter space. Therefore, one can start to constrain the seesaw 
parameters only after the discovery and subsequent determination of the SUSY parameters. 
In fig. [4] we show a similar graph but for the two generation model. The interesting point is 
that despite the fewer parameters one still has sufficient freedom to suppress the rare lepton 
decays. One the one hand, this shows the need to determine not only the differences of the 
neutrino masses squared but also the absolute neutrino mass scale or in other words the 
mass of the lightest neutrino, as a non-zero value of the latter would rule out the minimal 
two generation model or requires its extension by non-renormalizable operators. On the 
other hand it implies that for the exploration of the LHC phenomenology it is sufficient to 
study the simpler two generation model. 
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FIG. 5. BR(/i — > erf) (blue), BR(r — > e~f) (red) and BR(t — > jjL"f) (yellow) as function of sin<^> for 
two 24-plets with masses M\ = 5 • 10 14 GeV and M2 = 5 • 10 13 GeV on the left panel (right panel 
vice versa Mi > M2), 5 = and the SUSY parameters as in eq. (23). 



Up to now we have assumed that the .R-matrix is the unit matrix. In fig. [5] we study 
the dependence on the i?-matrix in the two-generation model. We fix the 24-plet masses 
to 5 • 10 13 GeV and 5 • 10 14 GeV and vary sin(0). Note that in both cases we have taken 
cos0 > 0. Instead of taking Mi > M2 we could have taken cos0 < in the second plot. As 
expected, we find that variation of the i?-matrix provides additional possibilities to suppress 
BR(/i ej) below the current experimental bound. Moreover, our results show that one 
never can exclude this class of models by these measurements as with a sufficient tuning of 
parameters one can always avoid the bounds. 
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FIG. 6. Comparison of the two-body decays li — > ^7 (left panel) and the-three body decays /j — > 3lj 
(right panel) for variation of sin(^i) at normal neutrino mass hierarchy, Dirac phase 5 = and a 
24-plet hierarchy: M ± = 10 15 GeV, M 2 = 10 14 GeV and M 3 = 10 13 GeV; k, lj = n, e (blue); r, e 
(red); r, \x (yellow). 

For completeness we compare in fig.[6]the branching ratios of the two body decays li — > Ijj 
to the ones for the three body decays U — > 3lj in the three-generation model. Similar to 
the seesaw type I case j32] we see that both decay classes show the same dependence on 
the underlying parameters, since in case of the three body decays the photon contribution 
dominates. We have checked that this also holds for the two-generation model. 
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B. Testing flavour structures at the LHC 



We have seen in the previous section that one can choose the seesaw parameters such 
that the experimental constraints on the rare leptons decays are fulfilled. In this section we 
address the question if there are any possibilities to test these models for such parameter 
choices at the LHC. As we will demonstrate there are indeed favorable SUSY parameter 
regions where one can observe the corresponding flavour violating decays of supersymmetric 
particles. 

The branching ratios of the lepton flavour violating decays of sleptons and neutralinos are 
governed by the same entries in the slepton mass matrix as the rare lepton decays, i.e. the 
ones given in eqs. (14) and (15). Therefore both classes of decays show the same dependence 
on the seesaw parameters. 

At the LHC one has to study cascade decays containing sequences of the form x\ ~>* 
l^lj — > iflfxi with i 7^ j |3"o] I77H83| . Moreover, the nature of the neutralinos should be 
dominantly gaugino like and the mass difference should be small enough to suppress the 
decay into h°. This requires a certain hierarchy between the neutralino mass parameters 
and the slepton mass parameters which is roughly given by \fi\ M2 > mj > Mi where the 
ordering of M\ and M2 can be interchanged. 

As the scaling of the lepton flavour violating decays of SUSY particles is similar to the 
one of the rare lepton decays in this class of models we use the following strategy to enhance 

the rates for X2 ~^ z^t^Xii X2 ~* ^ ±r=F Xi an d X2 ~~ /^ ±r=F X? : f° r a given point in the SUSY 
parameter space we choose the seesaw parameters in the following way: we fix the i?-matrix 
to be either 11 or as in eq. (12), depending whether we work in the two- or three-generation 
seesaw model. Next we fix the relative size of various entries of Yw such that the neutrino 
mixing matrix is tribimaximal. Note, that a non-zero #13 changes the neutralino branching 
ratios only slightly and, thus, its effect can be neglected her. In the third step both Yw and 
M w are rescaled until the correct neutrino masses are obtained and 10 12 • BR(/i —> ej) is in 
the interval [2.2, 2.4]. In this way one obtains the maximum rate for the decay X2 ~^ ^ ± / i=F Xi 
which is the cleanest at the LHC [80, 81j. With a further variation of the entries in Yw one 
could increase the final states containing a r lepton. However, we have checked for couple 
of points in the SUSY parameter space that this would only lead to an relative increase of 
about ten per-cent for the corresponding rates. We have not pursued this further as this 
is at most of the order of the expected theoretical uncertainty on the SUSY cross section, 
which is about 10-20 per-cent, see e.g. [84T487] and refs. therein. In the following examples 
we have checked that the bounds on SUSY particles are fulfilled [88 -91 J . 

The branching ratios of the lepton flavour violating decays can reach up to a few per- 
cent as shown in fig. [7j The structure of the RGEs implies that the three heaviest sleptons 
are essentially II even though there can be sizeable mixing between the stau states. The 
latter mixing is the main source of the LFV decays for Mi/ 2 < 1550 GeV where only the 
three lightest sleptons appear in the x° decays. The hierarchy BR(%2 Xi T f i ) > BR(%2 
XiA* e ) > BR(x° — > Xi Te ) is a consequence of the structure of Y w needed to explain the 
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The seesaw parameters are fixed as explained in the text. 



neutrino data. The change of the spectrum has two main sources: (i) M\/2 enters the RGEs 
for the slepton mass parameters, (ii) The requirement that BR(/i — > cy) to be in above 
interval implies that the seesaw scale becomes a function of Mu2- Changing the seesaw 
scale has a major impact on spectrum as discussed in detail in ref. Similar features 

show up in the 2-generation model as exemplified in fig. [8] where we show the LFV X2 decay 
branching ratios as a function of Aq. In this model one can find LFV branching ratios of up 
to 10 per-cent. The main reason for this are the different kinematics for the same mSUGRA 
input because changing the number of seesaw particles implies changes in the RGEs of the 
slepton and gaugino mass parameters as discussed above. 

We concentrate in the following on the two generation model as here the signal is some- 
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FIG. 9. Grid calculation of a(xf) X BR(x° ->■ X?A*e) (blue), a(xjj) x BR(x^ -»■ x?re) (red) and 
^Ow) x BR(x2 — ^ X? 7 "/ 2 ) (yellow) in femtobarn over M1/2 for different values of mo; ^4o = 0, 
tan/3 = 10 and fi > 0. The seesaw parameters are fixed as explained in the text. 
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FIG. 10. Invariant mass distributions for the signal and SUSY background for final states containing 
fj,, e (left plot) or r, // (right plot), missing transverse energy and at least two jets in the final state 
for a luminosity of 300 fb -1 , m = 50 GeV, M 1/2 = 1484 GeV, A = 0, tan/3 = 10 and \i > 0. 



what larger than in the three generation model. At the LHC X2 * s mainly produced in the cas- 
cade decays of squarks and gluinos. In fig. [9] we show ax BR as a function of Myz fixing the 
other parameters for two values of mo, Aq = 0, tan /3 = 10, fi > and ^/s = 14 TeV. In ad- 
dition, we used M X = M 2 = 2.5- 10 13 GeV as well as lV,n = ^V,i2 = -Y w ,ia = -5.252- 10" 2 , 
Yw,2i — and 1V,22 — Yw,23 — —1-547 ■ 10 _1 . Here we have summed over all possibilities to 
produce squarks and gluinos and we require that the two leptons from X2 are the only ones 
in the event. For the calculation of the cross section we have used the LHC-FASER package 
[92| 193]. One sees that the signal cross section before putting any cuts can be at most a few 
fb which gives at most a few hundred events even for an integrated luminosity of 300 fb _1 . 

This naturally leads to the question if such a signal can be observed at all. For this 
reason we have performed a Monte Carlo study at the parton level taking m = 50 GeV, 
M1/2 — 1484 GeV, A = 0, tan/3 = 10 and fi > corresponding to maximum of the 
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signal in fig. [9] For the generation of the events we use WHIZARD |71j. The corresponding 
signal cross sections are 1.4, 0.8 and 3.8 fb for the final states containing e/i, er and fir, 
respectively. Related studies have been performed in refs. [77, 80], [94] where it has been 
shown that one can reduce the SM background sufficiently. However, to our knowledge the 
SUSY background has not yet been taken into account. This will be also considered here. 
The main SM background are due tt, VV and Vjj (V = W, Z) production. In ref. [95J, 
where a detector study for the \xr channel has been performed for yfs = 10 TeV, it has been 
shown that the SM background can reduced significantly by requiring a cut on the missing 
transverse energy Ifix > 140 GeV and a cut on the effective mass M e ff > 400 GeV where 



The first sum is over the transverse momentum of the four hardest jets and the second one 
over the transverse momentum of all leptons. We have adjusted these cuts for the case 
a/s = 14 TeV and use IpT > 150 M e ff > 1200 GeV. Moreover we require that the event 
contains exactly two leptons and no 6-jets. This reduces the SM-background to a negligible 
level. The main SUSY background is due to charginos and W^-bosons produced in the SUSY 
cascade decays. In contrast to the signal these events stem in general from cascade decays 
of different squarks and/or gluinos. Therefore, if one plots the differential cross section as 
a function of the invariant lepton mass mu> = a/ (pi + Pi') 2 one gets a triangle for the peak 
and a flat distribution from the background. We have simulated the combination of signal 
with SUSY background using the dominant production mechanism which is in this case 
squark-squark production as the squarks are much lighter than the gluino yielding about 
80 per-cent of the total cross section. The results for an integrated luminosity of 300 fb -1 



are shown in fig. 10 where we have cut the range of mw at 500 GeV even though the SUSY 
background continues flat until about 1 TeV. As can be seen, one gets approximately the 
triangular shape of the signal with the edge at 

< = — * — i -y — ^ (24) 

where the lepton masses have been neglected. The edge clearly indicates the consecutive two 
body decays giving a first hint on the mass ordering. As the sleptons have different masses, 
they give somewhat different values for the edges which are collected in Tab. [T| Figure 10 



clearly shows that in this case the SUSY background is negligible compared to the signal. 
Note that the light sleptons hardly contribute to the signal as argued above and, thus, the 
edges are essentially due to the heavier sleptons. In case of the t/j, final state the two edges 
could be guessed but it will require high luminosity and a finer binning to disentangle the 
resulting double edge structure due to the contributions of the different sleptons [79] . 
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TABLE I. Edges, as given in eq. (24), of the invariant lepton masses due to the individual sleptons. 



The two neutralino masses are 344.6 GeV and 647.0 GeV 



slepton 


mass [GeV] 


m lv [GeV] 


h 


377.9 


213.0 




386.0 


233.9 


h 


621.9 


148.3 


h 


625.1 


139.3 


k 


625.9 


136.7 



IV. CONCLUSIONS 



We have studied supersymmetric variants of the seesaw type III model. At the electroweak 
scale the particle content is the same as in the MSSM. At the seesaw scale(s) the particles 
have been included in a 24-plet to ensure unification of the gauge couplings. In this way one 
ends up with a combination of seesaw type III combined with type I where the latter gives 
sub-dominant contribution if Si7(5)-GUT conditions for the corresponding Yukawa couplings 
are assumed. We have considered two variants of this model using either two or three 
generation of 24-plets. The latter case is heavily constrained by the experimental bound on 
fi — > cy. However, as we have shown there are various ways to obtain cancellations between 
different contributions so the bound can be respected: here the Dirac phase of the neutrino 
sector enters as well as the mass hierarchy of the seesaw particles and their mixing properties. 
Even though the measurement of the reactor angle #13 gives an additional constraint, the 
model still has sufficient many parameters to be consistent with all experimental data. In 
the two generation model the constraints due to the rare lepton decays are less severe and 
can be more easily accommodated. 

We have also investigated the question to which extent lepton flavour violating signals 
can be seen at the LHC The current experimental bounds on SUSY particles imply that 
within a unified model such as the mSugra squarks and gluinos must be in the TeV range. 
As the main signal is in the cascade decays of these particles one gets at most a few fb for 
the signal. The corresponding part of the parameter space is for small m and large rrii/2 if 
the seesaw parameters are chosen such that BR(/i — > cy) is close to its experimental bound. 
One can turn this around: if the bound BR(yU — > e'j) is increased by an order of magnitude 
than it is rather unlikely that LHC finds LFV in SUSY decays in this class of models. 
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